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We evaluate the decay branching ratios of xa — >■ PS, in a quark model parametrization scheme, 
where P and S stand for pseudoscalar and scalar meson, respectively. An interesting feature of this 
decay process is that the cc annihilate via the pQCD hair-pin diagram is supposed to be dominant. 
^— «j , Hence, this decay process should be sensitive to the quark components of the final-state light mesons, 

■ and would provide a great opportunity for testing the mixing relations among the scalar mesons, 
i.e. /o(1370), /o(1500) and /o(1710), by tagging the final state pseudoscalar mesons. 
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CS) ■ I. INTRODUCTION 

■ The usefulness of charmonium hadronic decays into light mesons is that this transition occur via a gluon-rich 
Oh process. The initial charm and anti-charm will annihilate into gluons and light quarks will be produced in the 
^ ■ final state through the hadronization of intermediate gluons. For the interest of studying the structure of the final 
Qj _ state light mesons, especially to probe possible existence of glueball states, the hadronic decays of charmonium system 

. provides an ideal platform on which the production of exotic states can be correlated with a relatively well- understood 
state. Such a tag sometimes exposes unexpected phenomena for which various possible mechanisms can be examined. 
I ! Review of heavy quarkonium dynamics can be found in Refs. [E Q. 
^ ■ In the charm sector there are several recent observations pertaining to scalar mesons in charmonium decays which 
OO ; turn out to be unexpected. One is the BES results for J/^p (jjf^ Q and fl, {i = 1, 2, 3 labels /o(1370), '/o(1500) 
OO ■ and /o(1710)), which show that the branching ratio of J/t/j — >■ (/)/o(1710) is smaller than that of J/ip — a;/o(1710). 
' In contrast, the branching ratio of J/ip a;/o(1370) is smaller than (/)/o(1370). The paradox arising here is that 
II /o(1710) dominantly decays into KK, hence is usually believed to have a large ss component. Therefore, one would 
PsJ naturally expect that the production of /o(1710) recoiled by 4> meson should be favored than by the oj due to the 
I OZI rule. Similar paradox occurs with the /o(1370) which is strongly coupled to iir and believed to be dominated by 
non-strange qq component. Efforts have been made in the literature to explore the scalar meson mixings As a 

" . I [ result, the light quark components in scalar mesons should be studied seriously in order to understand the properties 
J> ■ of these scalar mesons. 

The other experimental observation comes from CLEO-c and BES measurements p^ - flTj of Xc0,2 decays into 
^ . meson pairs. In Xco.2 ^ SS it was found by BES [13] that one of the largest branching ratios is from Xco ~^ 
' /o(1370)/o(1710). Again, the paradox arises from the OZI rule expectation. Compared with the decay channels of 
Xc0,2 — >■ PP and VV, one can see the better respect of the OZI rule [l^, [l^, where the DOZI processes are much 
suppressed. As a consequence, the leading transition amplitude is given by the singly OZI disconnected (SOZI) 
transitions. Therefore, the decay branching ratios for Xc0,2 ^ PP and VV still fit in the pattern of SU(3) flavor 
symmetry. In contrast, in the channel where scalar isoscalar /g states are produced. There are obvious deviations 
from the expectation of the OZI rule. As pointed out in Ref. [ISj], such a deviation could be a strong evidence for the 
glueball-gg mixing in the scalar isoscalar wavefunctions. Nevertheless, the scalar production is correlated with the 
large OZI-rule breaking in the DOZI processes as a common dynamic feature. 

Compared to Xc0,2 decays into two gluons, the cc annihilations in Xci decays would be dominated by the pQCD 
hair-pin diagram as shown in Fig.[TJa). The connected diagram (Fig.lTfb)), which generally plays a dominant role in 
XcQ.2 — ^ 'light hadrons', turns out to be suppressed by the Landau- Yang theorem [20| since the favored kinematics in 
Fig. mb) would be via two on-shell gluons. Empirically, the Landau- Yang suppression would occur in Figs. [2c)-(e), 
hereafter we refer to it as gluon loop diagrams. The following features should be recognized: i) The contributions from 
the DOZI diagram, i.e. Fig.[ljc), will be highly suppressed due to the large recoil momentum carried by the exchanged 
gluon between the quarks in the final state, ii) Figures [Ijd) and (e) may gain an enhancement by the gluon powers 
considering that the gluon couplings to the glueball generally do not pay a price. However, the suppression from the 
Landau- Yang theorem would be more essential since the kinematics of Figs.[T]^d) and (e), i.e. the loop diagrams, favor 
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to have those two gluons on sheU. Such a suppression, or in other words the dominance of the hair-pin diagram, is 
reflected by the total decay width of Xci, which is one order of magnitude smaller than that of the XcO- One notices 
that the Xc2 total width is compatible with that of the Xci- It can be explained by the Clebsch-Gordan coefficients 
when two vector gluons are coupled into a 2++ state. 

The three gluon annihilations of Xci into PS are also highly suppressed. One reason is because of the increase of 
the gluon powers. The other reason is that at least one of the final state meson will be produced via higher-twist 
components in the wavefunction as illustrated in Fig. [2l Based on the above considerations, it should be a reasonable 
approximation to treat the hair-pin diagram as the leading contribution to Xci ^ PS. 

In this work, we will show that by tagging the quark components of r/ and r/ in the final state, it is possible to 
investigate the quark components of the scalar mesons in xa — >■ PS. In particular, this process turns out to be 
sensitive to the glueball and qq mixing pattern. Therefore, it can be selective for different mixing schemes and serve 
as an alternative way to study the structure of scalar mesons. 

As follows, we first give details of the parametrization scheme in Sec. II. The numerical results are presented in 
Sec. Ill, and a brief summary is given in the last section. 



II. PARAMETRIZATION SCHEME 

A systematic parametrization scheme has been exploited for various charmonium hadronic decays E^, Uli j where 
the SOZI and DOZI processes can be parameterized out based on gluon counting rules. Those parameters can then be 
determined by experimental data from which predictions can be made for unmeasured channels. Early works based 
on similar parametrization can be found in Refs. [^-Q. In the case of Xci — ^ PS, as discussed in the Introduction, 
the dominant contribution is from Fig. [Ija) while the others are supposed to be strongly suppressed because of the 
Landau- Yang theorem. It makes the parametrization rather simple as we will detail below. 

For the production of isoscalar pseudoscalar mesons, since it has been well established that the glueball compo- 
nents inside 77 and rj' are rather small, it should be a good approximation to neglect their possible internal glueball 
components, and simply take the quark mixings in the SU(3) fiavor basis: 

77 = cos ap\nn) — siiiap\ss) , (1) 
7]' ~ aiTiap\nn) + cosap\ss) , (2) 

where ap = 9p + arctan\/2 and 0p is usually considered to be —22° ^ —13°. Here we adopt the commonly used 
value 6p = —19° to evaluated the branching ratios. With this mixing scheme, we eventually use the pseudoscalar 
mesons rj and 77' as the ffavor tag for the production of the quark components of the scalar mesons via Fig. [IJa) , while 
contributions from other processes can be neglected. 
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FIG. 2: The schematic Feynman diagrams of xa — >■ PS via three-gluon annihilations. 



The rich spectrum of scalar mesons in the mass region of 1^ 2 GeV has initiated a lot of studies of the scalar 
mesons including the search for the scalar glueball candidate (see Refs. [12, HI] and references therein). A broadly 
discussed scenario is the glueball-gg mixing among those three scalars, i.e. /o(1370), /o(1500) and /o(1710). In the 
SU(3) flavor basis, these states as the eigenstates of nn, ss and G components can be generally expressed as 

|/o(1710))\ / |G) \ /xi yi zi 

|/o(1500)) \^S\ \s7s) \ = \x2y2Z2\\ \ss) I , (3) 




|/o(1370)) J \ \nn) J \xs ys 



where Xi, yi and Zi are the mixing matrix elements determined by mixing mechanisms. Different models have 
different solutions for the mixing matrix critical difference among those mixing schemes focuses on 

the magnitude of glueball components inside /o(1710) and /o(1500), and in contrast, all those mixing schemes seem 
to agree that the /o(1370) is dominated by the nn component. Unfortunately, there still lack unique criteria for 
identifying the scalar glueball state and distinguish those mixing schemes. 

With Va standing for the potential of the SOZI process, a basic transition parameter g can be defined as 

g = {iqq)p{qq)s \ K | Xcl), (4) 

where q (q) is a non-strange quark (antiquark). Considering the SU(3) flavor symmetry breaking, which distinguishes 
an s quark pair production from the u/d quarks in the hadronizations, we introduce the SU(3) flavor symmetry 
breaking parameter R, 

^ ^ {{sq)p{qs)s I Vg I Xcl) j.^^ 
~ {{qq)p{qq)s I K | xci)' 

where i? = 1 is in the SU(3) flavor symmetry limit, while deviations from unity implies the SU(3) flavor symmetry 
breaking. In general, the value of parameter R is around R ~ ]k = 0.838, which provides a guidance for the 
SU(3) flavor symmetry breaking effects. For the creation of two pairs of ss via the SOZI process, the recognition of 
the SU(3) flavor symmetry breaking in the transition is 

^2 ^ {{ss)p{ss)s I Vg I Xcl) 

{{qq)p{qq)s I K | Xd) ' 

Following the above parametrization rule, we can write down the transition amplitudes for / = pseudoscalar {rj 
or f]') and scalar meson (/q, with i = 1,2,3 for /o(1710), /o(1500) and /o(1370), respectively) pair production as the 
following 

{vfo\Vg\Xci) = g{ziCosap - yiSmapR'^)T{P) (7) 
Wfa\yg\Xci) = giz^amap + y,cosapR^)J^iP) . (8) 
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For Xci decays into other channels with / 7^ 0, e.g. ttoo and KKq + c.c, the transition amphtudes are similar due to 
the exclusive contribution from Fig. [IJa); 

{^^%\Va\Xcl) = (^-^l^alXcl) - (^°«[!|K|Xcl) = 5-^(P) , (9) 

{K+K*-\V\xci) = {K-K*+\V\xci) - (X°Xo*VlXci) = K*°\V\xci) = gR^P) ■ (10) 

In the above equations, -^(P) is a commonly used form factor defined as follows, 

J-2(P) ^ |P|2'cxp(-PV8/32), (11) 

where P and I are the three-vector momentum and the relative orbit angular momentum of the final-state mesons, 
respectively, in the Xci rest frame. We adopt (3 — 0.5 GeV, which is commonly adopted in the literature IsL's']. At 
leading order the decays of Xci — ^ PS are via P-wave, i.e. I ~ 1. This form factor accounts for the size effects arising 
from the spatial wavefunction of the initial- and final-state mesons in the hadronizations. 



III. NUMERICAL RESULTS 

So far, the only available experimental information was given by BES-II [1^ HBl: i-e. BR{xci ^ i^}°(1430)XO -I- 
c.c. ^ K'^K+TT- + c.c.) < 8 X 10"'' and BR{xci ^ Kj+{U30)K~' + c.c. K^K+tt- + c.c.) < 2.3 x lO'^, where 
the statistics were limited and the spin of ii'}*'(1430) has not been determined. In this process both i4'Q(1430) and 
isr|(1430) may have contributions. However, it can be argued that the production of if|(1430)A' -f c.c. is relatively 
suppressed by the Clebsch-Gordan coefficients when those final-state fields are coupled together via a P wave. By 
assuming that the measured branching ratio is dominated by /Cq (1430)A'-|-c.c., it will set up upper limits of branching 
ratios of other decay channels. In another word, we can normalize the branching ratios of other decay channels to 
the i^Tp" (1430)7^° -I- c.c. channel (we take the lower limit as a conservative estimate), and inspect the variation 
of branching ratio fractions within a range of the SU(3) flavor symmetry breaking parameter R. Notice that the 
/•s:o*(1430) decay is nearly saturated by the Kt: channel with Pi?(X^(1430) Kir) = (93± 10)%, which is larger than 
BR{K2{li30) 3> Ktt) = (49.9 ± 1.2)% ^2^. We can simply take into account the charge conjugate by a factor of 3/2 
for K*°{U30)K" + c.c. -> K°K+Tr- + c.c. to obtain, BR{xci ^ ii'}°(1430)/-i:" + c.c.) < 1.2 x lO'^ as the upper limit. 

By studying this observable for different scalar meson mixing schemes, we can identify criteria for the determination 
of the quark contents of those scalar mesons. In turn, we expect to gain insights into their structures. 

From Eqs. ([7]) and ([T0|. we can take the branching ratio fraction 

BRixci-^vfa) ^ |P,) I (zi cos gp-y, sin apP^)^J'^(P^) 

- BRixci ^ KOR*" + c.c.) 2\Pk\R^T^Pk) ' ^ ' 

where the factor 2 in the denominator is due to the charge conjugate factor for K'^Kq'^ + c.c. final states. Similarly, 
for the rj' recoiling /q we can define 

BR{xci ^ v'f^) _ IPy I {Zj sin ap + cos apR^fJ^^ [P^, ) 

- Bi?(xci ^ KOK*° + c.c.) 2\Fk\R'T^Pk) ' ^ ' 



and for the tt'^Oq channel we have 



BRixci^TT^a!?,) _ \P^\T^iP^) 
- BRixci ^ K^K*" + c.c.) 2\Pk\R^J^HPk) ' ^ ^ 



For Xci K^Kq^ + c.c, the only parameter is R for which R ~ jirl !k ~ 0.838 is commonly adopted. Applying 
the experimental upper hmit, BR{xc\ i^}(1430)°i^" + c.c.) < 1.2 x 10~^ [H, [2^, we can determine the basic 
transition strength g = 2.56 x 10~^ via 

r(Xci ^ ifj(1430) + c.c.) = yI-kAP ' ^ ^ 

Also, the branching ratio of Xci ^^Q^o can be estimated via Eq. p^ . i.e. BR^^^^^OgO^ < 8.64 x 10~^. 

Now we focus on Eqs. ((T^ and ((T5| to extract information about the scalar meson structures. The ratios 7^ and 7^' 
are now explicit functions of the qq mixing elements. We will analyze three typical mixing schemes in the literature. 
Predictions of the ratios would set up criteria for future experimental tests of those mixing scenarios. 
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FIG. 3: The branching ratio fractions BR{xci viv') fo) / B R(xci K}{l'iSO)'^ K'^ + c.c.) as a function of the parameter R 
is presented in this figure. In this case, the mixing matrix is used for scalar meson mixing. The solid line, dashed line and 
dotted line are the branching ratio fractions of r;/o(1710), r;/o(1500) and r;/o(1370) in Diagram (a). The cases in Diagram (b) 
are for r77o(1710), r;7o(1500) and t?7o(1370). 



Scheme-I: 

A systematic study by Close et al based on a perturbation transition mechanism [SiS] determines the glueball-gg 
mixing matrix Q, 

/ 0.36 0.93 0.09 \ 
Si = -0.84 0.35 -0.41 . (16) 
\ 0.40 -0.07 -0.91 / 

In this scheme the /o(1710) is dominated by the ss component, but its glueball component is also sizeable. In 
contrast, the /o(1500) is dominated by G with a sizeable ss component. The /o(1370) is found dominated by the 
nh (= {uu + dd)/^/2). 

Consider a small variation of the parameter R, we plot the branching ratio fractions 7^ and 7^' in Fig. [3] with the 
mixing matrix elements from Eq. (I16p for the mixing Scheme-I. It is interesting to learn the following points about 
this scenario: 

• Comparing the branching ratio fractions 7,, and 7,,' , we notice that the relative phases between the mixing 
matrix elements critically determine the production strengths of those scalars when recoil 77 or 77'. 

• Since the nn dominates the /o(1370) wavefunction and the ss is negligibly small, the production rates of /o(1370) 
are predicted to be larger than other channels in the vicinity of i? ~ 0.838. 

• For the /o(1710) its production in association with 77 is relatively suppressed by the significant cancelation 
between the ss and nn components. In contrast, its production with rj' is much more enhanced. 

• The ss and nn components are compatible in the /o(1500), but out of phase. As a consequence, the production 
of the /o(1500) seems to be unfavored in Xci PS for the mixing Scheme-I. 

Scheme-II: _ 

The second mixing scheme was provided by Cheng, Chua and Liu 10] based on the lattice QCD (LQCD) quenched 
calculations by Lee and Weingarten [9(1 . It was found that /o(1710) is dominated by the glueball component, while 
/o(1500) and /o(1370) are dominated by ss and nn, respectively. This scenario is different from Ref. [I3] where the 
glueball-dominant state is the /o(1500). The mixing scheme of Ref. [l^ is similar to that of Rcf. Q. Therefore, we 
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FIG. 4: The notations are similar to Fig[3]with the mixing matrix 52. 
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FIG. 5: The notations are similar to Fig[3]with the mixing matrix S^a- 



do not show the numerical survey of these two solutions, but only present the results for the typical solution from 
Ref. fi 



0.859 0.302 0.413 
S2 = \ -0.128 0.908 -0.399 
-0.495 0.290 0.819 



(17) 



The predicted branching ratio fractions 7^ and 7^' in terms of R are plotted in Fig. |4l In comparison with Fig. |3l 
the predicted decay pattern is obviously different. In particular, one notices the strong suppression of 77/0 (1710) in 
comparison with 77/o(1500) and 77/o(1370) in the vicinity of i? ~ 0.838. In the 77' production channels, the 77'/o(1710) 
and 77'/o(1500) decays are relatively suppressed in comparison with 77'/o(1370). 
Scheme-Ill: 
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FIG. 6: The notations are similar to Fig|3]with the mixing matrix Ssb- 



Giacosa et al. obtained four possible solutions by fitting the masses and decay widths of those three /o states in 
an effective chiral approach [24l |. Their two typical solutions, i.e. S^a and S^t were obtained without direct glueball 
decays, and have the following expressions. 



and 



Sab — 



-0.06 0.97 -0.24 
0.89 -0.06 -0.45 
0.45 0.24 0.86 



-0.68 0.67 -0.30 
0.49 0.72 -0.49 
0.54 0.19 0.81 



(18) 



(19) 



These two solutions were extracted with the OZI-rule violation parameters r — 1.93 ± 0.29 and —2.07 ± 0.79 
determined in J/ip — > ipf^ and cj/q, and they are both of order of one. The main difference between these two 
solutions lies in the different prescriptions of the glueball and ss contents for the /o(1710) and /o(1500), respectively. 
Compared with mixing scheme-I and II, it shows that S^a is similar to of Ref. Q, but the second solution S^b is 
quite different. 

In Figs. [S] and [HI the predicted branching ratio fractions are presented for those two solutions S3 (I and 5*3^, respec- 
tively. Their different prescriptions lead to drastic changes of the production rate for the /o(1500) around R — 0.838. 
In the mixing scheme of S^a, the decay channel 77/o(1500) is suppressed in comparison with ?7/o(1710) and 77/o(1370), 
while in 5*36 it is strongly enhanced to be larger than the other two channels. The R dependence of 7,, and 7,,' can 
also be observed in these two mixing schemes. 

Interestingly, in the vicinity of i? ~ 0.838 the decay patterns illustrated by those mixing schemes (four different 
mixing matrices) can still be distinguished. As mentioned earlier that the structure of mixing matrices, S^a and 6*1, 
are similar to each other, we can see that their predictions for the production of /o(1710) are quite different. For 
instance, it shows that the production of /o(1710) in association with r/ is more favored than with 77' in S^a, and it is 
opposite in 5*1. 

In TablelH we list the branching ratios of Xci v{v')fo with R — 0.838 as a predictions from those mixing schemes. 
Combining what illustrated in Figs. [5]|ni we learn the following points concerning the production of the scalars: 

• The /o(1370) is dominated by the nn component in all those schemes. Nevertheless, all those mixing schemes 
find that the relatively small ss component is in phase to the dominant nn. As a consequence, the predicted 
branching ratios of both ?7/o(1370) and r//o(1370) turn out to have a stable behavior. 
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TABLE I: The upper limits of the branching ratios Xci rj fo is presented with the SU(3) breaking parameter R = /,r//i^ 
0.838. The experimental value BR{xci -i> /C}(1430)"^" + c.c.) < 1.2 x 10~^ 26] is used to predict the upper limits. 



BRjxci PS) Si S2 Ssa Sat 

77/0(1710) 7.2 X 10"^ 3.3 X 10"^ 2.63 x 10"" 2.03 x 10~* 

r?/o(1500) 1.88 x 10"'' 3.23 x 10"* 9.69 x 10"^ 3.98 x 10"* 

r?/o(1370) 4.32 x 10"* 2.55 x 10"* 3.08 x 10"* 2.90 x 10"* 

?77o(1710) 2.10 X 10"* 1.04 x 10"* 1.04 x 10"* 2.70 x 10"^ 

V/o(1500) 4.89 X lO"'^ 8.27 x 10"^ 6.38 x 10"^ 1.09 x 10"'' 

r77o(1370) 2.51 x 10"* 3.21 x 10"* 3.17 x 10"* 2.61 x 10"* 



• For the /o(1500), sensitivities of the branching ratio fractions 7,, and 7,,' to its quark contents can be seen. The 
predicted branching ratios of Xci ~> '7/o(1500) are at the order of lO"'* which could be accessible at BES-III. 
One also notices that the 77'/o(1500) channel is relatively suppressed in all those mixing schemes since the nfi 
and ss have opposite signs in the mixing schemes 5*1, S2 and 5*36, while in S^a the nn and ss both are small. 

• The branching ratios of the ?7/o(1710) and 77'/o(1710) are also sensitive to the mixing schemes. Combine together 
the decay patterns of other scalars, it is possible to distinguish those mixing schemes in experiment. 

It should also be pointed that experimental analysis may become much more complicated due to the background 
contributions to the final states. For the PS channel, the final-state particles involve titttt, ijKK, ti'tttt, and rj'KK etc. 
As shown by the recent measurement from CLEO-c [131 that the rjTnr and ti'tttt channels may not be suitable for the 
search for the scalar meson signals due to large background contributions. In contrast, the t]KK and rj'KK channels 
may be more sensitive to the scalar meson productions. With the high statistics measurement at BEPC-II/BES-III, 
we anticipate that progress can be made in the study of the scalar meson spectrum. 

One should also be cautioned that the dominance of the hair-pin diagram is a crucial assumption in this study. A 
better respect of this assumption may be achieved in the bottomonium sector, namely, in Xbi — > PS. Unfortunately, 
there are no experimental data available at this moment for this channel. The future LHCb experiment may be able 
to provide additional information about the nature of those scalar mesons. 



IV. SUMMARY 



To summarize, we show that the decay of Xci — ^ PS could be an ideal channel for probing the quark contents 
of the those scalar mesons in the mass region of 1^2 GeV, i.e. /o(1370), /o(1500) and /o(1710). Because of the 
suppression from the Landau- Yang theorem to those gluon loop diagrams, this decay channel at leading order should 
be dominated by the pQCD hair-pin transition process. It thus allows us to tag the quark contents of the final state 
scalars by the quark components of the recoiled rj and rj' . A prediction for the upper limits of the branching ratios of 
Xci vfo ^-i^d Ty'/o can be made with the available experimental data and based on different scalar mixing schemes 
in the literature. It can be expected that a precise measurement of Xci — > PS will be able to distinguish those model 
prescriptions for the glueball-gg mixing scenario and provide further evidence for the scalar glueball candidates. 
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